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Abstract 



The post-Newtonian gravitoelectric secular rate of the mean anomaly 
M. is worked out for a two-body system in the framework of the Gen- 
eral Theory of Relativity. The possibility of using such an effect, which 
is different from the well known decrease of the orbital period due to 
gravitational wave emission, as a further post-Keplerian parameter in 
binary systems including one pulsar is examined. The resulting effect 
is almost three times larger than the periastron advance u>. E.g., for 
the recently discovered double pulsar system PSR J0737-3039 A+B 
it would amount to -47.79 deg yr _1 . This implies that it could be 
extracted from the linear part of a quadratic fit of the orbital phase 
because the uncertainties both in the linear drift due to the mean 
motion and in the quadratic shift due to the gravitational wave are 
smaller. The availability of such additional post-Keplerian parame- 
ter would be helpful in further constraining the General Theory of 
Relativity, especially for such systems in which some of the other post- 
Keplerian parameters can be measured with limited accuracy. More- 
over, also certain pulsar-white dwarf binary systems, characterized by 
circular orbits like PSR B1855+09 and a limited number of measured 
post-Keplerian parameters, could be used for constraining competing 
theories of gravity. 



1 The post-Newtonian rate of the mean anomaly 

According to the Einsteinian General Theory of Relativity (GTR), the post- 
Newtonian gravitoelectric two-body acceleration of order 0(c -2 ) (1PN) is, in 
the post-Newtonian centre of mass frame (see [3] and, e.g., |12j and references 
therein) 
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where r and v are the relative position and velocity vectors, respectively, G 
is the Newtonian constant of gravitation, c is the speed of light, mi and 777-2 
are the rest masses of the two bodies, m = mi + 777,2 and v = 77717772/m 2 < 1. 
The orbital phase can be characterized by the mean anomaly M defined 

as 

M = n{t-T ), (2) 
where the unperturbed mean motion 77 is defined as 

n = Tb , (3) 

Pb is the anomalistic period, i.e. the time elapsed between two consecutive 
pericentre crossings, which is 2iry/ a^/Gm for an unperturbed Keplerian el- 
lipse, and To is the date of a chosen pericentre passage. 

The variation of the mean anomaly can be written, in general, as 

dM (Pb\ ( . ^ 2vr fdT \ 

— = n-2^j(t-T )--^— j. (4) 

The second term of the right-hand side of (Jl|) accounts for any possible 
variation of the anomalistic period. The third term, induced by any small 
perturbing acceleration with respect to the Newtonian monopole, whether 
relativistic or not, is the change of the time of the pericentre passage, which 
we will define as 

dt P b \dt J ' U 
It can be calculated with the aid of the Gauss3 perturbative equation [13] 
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where R is the radial component of the disturbing acceleration, a, e, i, f2 and 
ijj are the semimajor axis, the eccentricity, the inclination, the longitude 
of the ascending node and the argument of pericentre, respectively, of the 
orbit. In order to obtain the secular effects, we must evaluate the right- 
hand-side of © on the unperturbed Keplerian ellipse and, then, average the 
result over one orbital revolution. In regard to the physical interpretation 
of d£/dt, the apsidal line can be crossed at a different time with respect to 



For a different approach based on a modified form of the Lagrange planetary equations 
see 12]. 
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the unperturbed Keplerian case either because there is an additional radial 
acceleration which accelerates/decelerates the moving particle or because 
the pericentre is no more fixed in space. The first case is accounted for by 
the first term of the right-hand side of ©, while the second term of the 
right-hand side of © explains the second case. Indeed, for, e.g., an inward 
radial acceleration which increases the gravitational strength the particle 
moves more rapidly and the pericentre, assumed fixed in space, is reached 
in advance with respect to the unperturbed case: dTo/dt < 0. If the angle 
u + f2 cos i, which defines an angular variable around the axis of the orbital 
angular momentum, increases, the passage at pericentre will occur later with 
respect to the unperturbed case: dTo/dt > 0. 

We will now consider (pQ) as perturbing acceleration. Let us start with 
the first term of the right-hand-side of ©. By defining 



A 



B 



C 



(Gmf 



(4 + 2i/), 



Gm, 



l + 3i/), 



(7) 



^(4-f), 



it is possible to obtain from ([I]) 
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Now the term —2Rr/na 2 , with R given by ([8]), must be evaluated on the 
unperturbed Keplerian ellipse characterized by 
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where / is the true anomaly, and averaged over one orbital period by means 
of 

dt r 2 df 



Thus, 
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In the expansion of r in (jlip the terms of order C(e 4 ) are retained. The 
final result is 

2 r\ nGm 
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The post-Newtonian gravitoelectric secular rate of pericentre is indepen- 
dent of v and is given by the well known formula 
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while there are no secular effects on the node. 

The final expression for the post-Newtonian secular rate of the mean 
anomaly is obtained by combining (jl2H - A14H and by considering that, for a 
two-body system, it is customarily to write 
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where M = w/Mq is the sum of the masses in units of solar mass and 
T = GM /c 3 = 4.925490947 x 10~ 6 s. It is 
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Note that (fT6|) is negative because (fT7|) is always positive; thus the crossing 
of the apsidal line occurs at a later time with respect to the Kepler-Newton 
case. 

Note that, for v — * 0, i.e. mi <C m,2, (I16j) does not vanish and, for small 
eccentricities, it becomes 

(18) 
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which could be used for planetary motion in the Solar System [7j. E.g., for 
Mercury it yields a secular effect of almost -130 arcsec cy" 1 . It is important 
to note that the validity of the present calculations has also been numerically 
checked by integrating over 200 years the Jet Propulsion Laboratory (JPL) 
equations of motion of all the planets of the Solar System with and without 
the gravitoelectric 1/c 2 terms in the dynamical force models [5] in order to 
single out just the post-Newtonian gravitoelectric effects. They fully agree 
with (|18p (E.M. Standish, private communication, 2004). Another analytical 
calculation of the post-Newtonian general relativistic gravitoelectric secular 
rate of the mean anomaly was performed [TJj in the framework of the La- 
grangian perturbative scheme for a central body of mass M-test particle 
system. The author of [13] starts from the space-time line element of the 
Schwarzschild metric written in terms of the Schwarzschild radial coordi- 
nate r . Instead, (H|) and the equations of motion adopted in the practical 
planetary data reduction at, e.g. JPL, are written in terms of the stan- 
dard isotropic radial coordinate r related to the Schwarzschild coordinate 
by r = r(l + GM /2c 2 r) 2 . As a consequence, the obtained exact expression 
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contrary to the pericentre case, agrees neither with (|18|) nor with the JPL 
numerical integrations yielding, e.g., a secular advance of +42 arcsec cy -1 
for Mercury. For a better understanding of such comparisons, let us note 
that both the numerical analysis by Standish and of the author of [H] are 
based on the = case; n gets canceled by construction in the Standish 
calculation, while in [14] the numbers are put just into (I19p . which is the 
focus of that work. 



1.1 Testing gravitational theories with binary pulsars 

In general, in the pulsar's timing data reduction process^] five Keplerian 
orbital parameters and a certain number of post-Keplerian parameters are 
determined with great accuracy in a phenomenological way, independently 
of any gravitational theory [191 [T6] . The Keplerian parameters are the pro- 
jected semimajor axis x = asinz/c, where i is the angle between the plane 
of the sky, which is normal to the line of sight and is assumed as reference 
plane, and the pulsar's orbital plane, the eccentricity e, the orbital period 
Pfe, the time of periastron passage To and the argument of periastron u>o at 

2 For all general aspects of the binary pulsar systems see [191116] and references therein. 
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the reference time To. The most commonly used post-Keplerian parameters 
are the periastron secular advance uj, the combined time dilation and gravi- 
tational redshift due to the pulsar's orbit 7, the variation of the anomalistic 
period Pi,, the range r and the shape s of the Shapiro delay. These post- 
Keplerian parameters are included in the timing models |19[ [T5] of the so 
called Roemer, Einstein and Shapiro Ar, Ag, As delays^] occurring in the 
binary pulsar systemEl 

Ar = x sin uj [cos E — e(l + 5 r )] + xcosusinE^l — e 2 (l + 5g) 2 , 

< Ae = 7sini?, 

As = — 2r ln{l — ecos E — s[sinw(cos E — e) + y/l — e 2 cos a; sinE 1 ]}, 

(20) 

where E is the eccentric anomaly defined as, E — e sin E = M . cos E and 
sinE 1 appearing in (|20p can be expressed in terms of Ai by means of the 
following elliptic expansions [T7j 

cos E 

< 

sini? 

v 

where Jj(y) are the Bessel functions defined as 

f'TT 

7T J j (y)= cos(j0 -y sin 6) d8. (22) 
Jo 

The relativistic secular advance of the mean anomaly (|16p can be accounted 
for in the pulsar timing modelling by means of (|2ip . 

In a given theory of gravity, the post-Keplerian parameters can be writ- 
ten in terms of the mass of the pulsar m p and of the companion m c . In 
general, m p and m c are unknown; this means that the measurement of only 
one post-Keplerian parameter, say, the periastron advance, cannot be con- 
sidered as a test of a given theory of gravity because one would not have a 
theoretically calculated value to be compared with the phenomenologically 
measured one. In GTR the previously quoted post-Keplerian parameters 

3 For the complete expression of the timing models including, e.g., also the delays 
occurring in the Solar System due to the solar gravity see |19l I16j . 

4 The aberration parameters S r and 5g are not, in general, separately measurable. 
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It is important to note that the relativistic expression of Pf, in f|23j) . should 
not be confused with £ 
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of (|16|) . Indeed, it refers to the shrinking of the 

orbit due to gravitational wave emission which vanishes in the limit v —* 0, 
contrary to (|16|) which expresses a different, independent phenomenon. The 
measurement of two post-Keplerian orbital parameters allows to determine 
m p and m c , assumed the validity of a given theory of gravit}0. Such values 
can, then, be inserted in the analytical expressions of the remaining post- 
Keplerian parameters. If the so obtained values are equal to the measured 
ones, or the curves for the 2 + N, with N > 1, measured post-Keplerian 
parameters in the m p — m c plane all intersect in a well determined (m p , m c ) 
point, the theory of gravity adopted is consistent. So, in order to use the 
pulsar binary systems as valuable tools for testing GTR the measurement 
of at least three post-Keplerian parameters is required. The number of 
post-Keplerian parameters which can effectively be determined depends on 
the characteristics of the particular binary system under consideration. For 
the pulsar-neutron star PSR B1913+16 system [6] the three post-Keplerian 
parameters cj, 7 and P\, were measured with great accuracy. For the pulsar- 
neutron star PSR B1534+12 system [15] the post-Keplerian parameters re- 
liably measured are cj,7,r and s. For the recently discovered pulsar-pulsar 
PSR J0737-3039 A+B [lj system the same four post-Keplerian parameters as 
for PSR B1534+12 are available plus Pb and a further constraint on m p /m c 
coming from the measurement of both the projected semimajor axes. On the 
contrary, in the pulsar-white dwarf binary systems, which are the majority 
of the binary systems with one pulsar and present almost circular orbits, it 



5 This would still not be a test of the GTR because the masses must be the same for 
all the theories of gravity, of course. 
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is often impossible to measure Co and 7. Up to now, only r and s have been 
measured, with a certain accuracy, in the PSR B1855+09 system [8], so that 
it is impossible to use its data for testing the GTR as previously outlined. 



1.2 The secular decrease of the mean anomaly and the bi- 
nary pulsars 

Let us investigate the magnitude of the mean anomaly precession in some 
systems including one or two pulsars. 

For PSR B1913+16 we have [IS] m p = 1.4414M , m c = 1.3867M ,e = 
0.6171338, P b = 0.322997448930 d. Then, v = 0.2499064, F = 1.04459537192 
and £ = -10.422159 deg yr" 1 . For PSR J0737-3039 A we have [TO] 

GE 

m p = 1.337M , m c = 1.25OM , e = 0.087779, P b = 0.102251561 d, so that 

v = 0.249721953643, F = 0.946329857430. Thus, £ = -47.79 deg yr" 1 . 

GE 

This implies that the ratio of the post-Newtonian gravitoelectric secular rate 
of the mean anomaly to the mean motion amounts to ~ 10~ 5 . Let us see 
if such post-Newtonian shift is detectable from quadratic fits of the orbital 
phases of the form M = a® + b^t + cot 2 . For PSR B1913+16 the quadratic 
advance due to the gravitational wave emission over thirty years amounts 
to (P b = -2.4184 x 10~ 12 ) 

AM = -ir \ I (t - T ) 2 = 0.5 deg, (24) 



with an uncertainty 5(AM) fixed to 0.0002 deg by 5P b = 0.0009 x 10~ 12 . 
The linear shift due to (fTOl) amounts to 

(t - T ) = -312.6647 deg (25) 

GE 

over the same time interval. The uncertainty irj! n amounts to 1 x 10 -9 deg 
yr" 1 due to 5P b = 4 x 10~ 12 d. For PSR J0737-3039 A the gravitational 
wave emission over three years (P b = —1.20 x 10~ 12 ) induces a quadratic 
shift of 0.008 deg, with an uncertainty 5(AM) fixed to 0.0005 deg by 5P b = 
0.08 x 10" 12 . The linear shift due to CE]) amounts to -143.3700 deg over the 
same time interval. The uncertainty in n amounts to 7 x 10~ 7 deg yr _1 due 



AM = £ 



The sum of the masses and the semimajor axis entering n are determined from timing 
data processing independently of M itself, e.g. from the periastron rate and the projected 
barycentric semimajor axis. 
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to 5Pb = 2 x 10 10 d. Thus, it should be possible to extract £ 



from the 

GE 

measured coefficient bo; both the corrupting bias due to the uncertainties in 
the quadratic signature and the errors in n would be negligible. 

Measuring £ as a further post-Keplerian parameter would be very 
GE 

useful in those scenarios in which some of the traditional post-Keplerian 
parameters are known with a modest precision or, for some reasons, cannot 
be considered entirely reliabk^]. E.g., in the double pulsar system PSR 
J0737-3039 A+B the parameters r and 7 are measured with a relatively 
low accuracy [11]. Moreover, there are also pulsar binary systems in which 
only the periastron rate has been measured in this case the knowledge 
of another post-Keplerian parameter would allow to determine the masses 
of the system, although it would not be possible to constraint alternative 
theories of gravity. 



2 Conclusions 

In this paper we have analytically derived for a two-body system in eccentric 

orbits the secular variation £ yielding the post-Newtonian general rela- 

GE 

tivistic gravitoelectric part of the precession of the mean anomaly not due 
to the variation of the orbital period. In the limit of small eccentricities and 
taking the mass of one of the two bodies negligible, our results have been 
compared to the outcome of a numerical integration of the post-Newtonian 
general relativistic gravitoelectric equations of motion of the planets of the 
Solar System performed by JPL: the agreement between the analytical and 
numerical calculation is complete. Subsequently, we have investigated the 
possibility of applying the obtained results to the binary systems in which 
one pulsar is present. In particular, it has been shown that the variation of 
the orbital period P^ due to gravitational wave emission and the effect 

gw 

derived by us are different ones. Indeed, the post-Newtonian gravitoelectric 
precession of the mean anomaly, which is always negative, is related to the 
secular increase of the time of pericentre passage and occurs even if the or- 
bital period does not change in time. A quadratic fit of the orbital phase of 

the pulsar would allow to measure £ because the biases due to the errors 

GE 

in the quadratic shift due to P& and in the linear shift due of the mean mo- 



7 The measured value of the derivative of the orbital period A is aliased by several 
external contributions which often limit the precision of the tests of competing theories of 
gravity based on this post-Keplerian parameter |16| . 
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tion n are smaller. The use of £ as a further post-Keplerian parameter 

GE 

would allow to improve and enhance the tests of post-Newtonian gravity 
especially for those systems in which only few post-Keplerian parameters 
can be reliably measured. 
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